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O Abstract - Spatially homogeneous cosmological models reduce to Hamiltonian systems in 
O ■ a low dimensional Minkowskian space moving on the total energy shell H = 0. Close to the 
q ' initial singularity some models (those of Bianchi type VIII and IX) can be reduced further, 
in a certain approximation, to a non-compact triangular billiard on a 2-dimensional space 



of constant negative curvature with a separately conserved positive kinetic energy This 
^ type of billiard has long been known as a prototype chaotic dynamical system. These 
(3jT)! facts are reviewed here together with some recent results on the energy level statistics of 
J> ! the quantized billiard and with direct explicit semi-classical solutions of the Hamiltonian 
^ ■ cosmological model to which the billiard is an approximation. In the case of Bianchi 
■ type IX models the latter solutions correspond to the special boundary conditions of a 
'no-boundary state' as proposed by Hartle and Hawking and of a 'wormhole' state. 



INTRODUCTION 

The subject of chaos in cosmology is really an old one and predates the intense modern interest in chaotic 
dynamical systems. It started with the classical work of Belinsky et al 1969 |1[ and of Misner 1969 
who found that the evolution of certain homogeneous cosmological models back into the past towards 
the 'Big Bang' shows an oscillatory behavior, some aspects of which in modern terminology one has to 
call chaotic. Useful references for these cosmological models are the books by Ryan 1972 || and Ryan 
and Shepley 1975 [Q. The chaotic aspects of their evolution have been examined in numerous papers 
starting with the review by Barrow 1982 |J. Khalatnikov et al 1985 || give a detailed discussion. A 
dynamical systems approach was pioneered by Bogoyavlensky and Novikov 1973 p] and elaborated by 
Ma and Wainwright 1989 ||. Among the more recent works on this subject those of Pullin 1990 
and Rugh 1992 |TI| may be mentioned. Most recently a Painleve analysis was applied to this dynamical 
system by Cantopoulos et al 1993 |TT[ who found, surprisingly, that the system satisfies the Painleve 
conditions which are necessary (but not sufficient) for a system to be integrable. 
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Early on, again before the explosion of the interest in chaos took place, Misner 1972 \Y1\ in a seminal 
paper and his student Chitre 1972 W3\ in his thesis introduced coordinates in which the most interesting 
homogeneous cosmologies, those of Bianchi types VIII and IX, sufficiently close to the singularity reduce 
to a special two-dimensional billiard on a homogeneous space of constant negative curvature. The study 
of this kind of billiards is an even more ancient subject in mathematics, dating back to the work of 
Hadamard 1898 [|TJ[], Artin 1924 and Hopf 1937 jn|. An extensive review is due to Balazs and 
Voros 1986 |T7| . 

Given a dynamical system of Hamiltonian form the temptation to quantize seems to be irresistable. 
Hence the field of quantum chaos, which is concerned with the quantization of classically chaotic systems 
with the aim of understanding features of chaos and quantum mechanics, hence also the field of quantum 
cosmology, which starts with the quantization of spatially homogeneous cosmological models with the 
aim to achieve an understanding of quantum features of cosmology For the Bianchi type IX cosmology 
this was first done by Misner 1972 in [|12|]. However, the present interest in quantum cosmology was 
mainly triggered by the proposal of a specific initial condition, 

the 'no-boundary-state', by Hawking 1982 ||18[| , and Hartle and Hawking 1983 Quantized versions 
of the homogeneous Bianchi type IX cosmological models were also investigated in this context, in 
semi-classical WKB-type approximations, e.g. in papers by Hawking and Luttrell 1984 [pOR , and Moss 
and Wright 1985 |2l[ who study consequences of the 'no-boundary'condition and by Del Campo and 



Vilenkin 1989 |22], and Graham and Szepfalusy 1990 [23| who instead pose the boundary condition of an 
outgoing wave. From the point of view of quantum chaos, the Bianchi IX model has also been studied 
further in papers by Furusawa 1986 [24|] and Berger 1989 |25] who study the dynamics of wave-packets 
and by Graham et al. 1991 and Csordas et al. 1991 [fjTj who study level statistics in Misner's billiard 
limit of the model. More recently further detailed studies of a class of billiards including the special one 
relevant to the Bianchi VIII and IX model have appeared in the quantum chaos literature (Bogomolny 
29|, Bolte et al. 1992 0, Schmit 1991 H Eisele and Mayer 1993 ||). 



et al. 1992 

Finally, as the most recent development, exact solutions of quantized Bianchi IX models were found 
independently by Moncrief and Ryan 1991 0], Graham 1991 |34j], and Bene and Graham 1993 [35|. It is 
now understood that the solution found in the first two references describe virtual quantum wormholes 
(D.Eath 1993 [j36 , Bene and Graham 1993 f35j), while a further exact solution by Bene and Graham 
1993 gives the 'no-boundary-state' of the model in a simple analytical form. 

The present paper has the purpose to review some of these results. In the next section we establish 
the necessary background on spatially homogeneous cosmological models. Then we proceed to consider 
Misner's approximate asymptotic reduction of the Bianchi VIII and Bianchi IX models to 2-dimensional 
billiards on a space of constant negative curvature, and discuss some results following from this reduction. 
Next, the Bianchi VIII and IX models are quantized, giving rise to their Schrodinger equation, called 
Wheeler DeWitt equation in the present context. We then discuss some of the properties of the billiard 
approximation of the quantized model. Finally, we turn to exact solutions of the Wheeler DeWitt 
equation and their physical interpretation. 



SPATIALLY HOMOGENEOUS COSMOLOGICAL MODELS 

We consider the class of space-times such that through every point there passes a space-like 3-manifold 
which is left invariant by the actions of a three-dimensional Lie group (see e.g. Ryan and Shepley 1975 
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0). Such space-times are called spatially homogeneous. The three-dimensional Lie groups have been 
classified by Bianchi into nine different types ('Bianchi types'). They can be distinguished by the Lie 
algebra of the group-generators, called 'Killing vectors' in general relativity. It is convenient to use an 
invariant (non-coordinate) basis and its dual of basis 1-forms uf on the invariant 3-manifold. 
The space-time metric then takes the general form 

ds 2 = -dt 2 + gij(t)uj l uj j (1) 

where dt is the time-like one-form orthogonal to the invariant 3-manifold and t is the standard cosmic 
time coordinate. The components <7y(t) of the 3-metric then depend only on t. The basis 1-forms satisfy 

dJ = -CtjitJ A u k (2) 

(we recall the usual definition A u k = ® u k — u) k (£) up) where the C % jk are the structure constants 
of the homogeneity group. Of special interest here will be the Bianchi type VIII with 

— C 23 — +C 32 = +C 2 31 — —C 2 13 — +C 3 12 = ~C 3 21 = 1 (3) 

and in particular the Bianchi type IX with 

C l jk = tijk- (4) 

In the latter case the homogeneity group is SO (3), the invariant 3-manifold is compact and topologically 
a 3-sphere. This case therefore includes the closed Friedmann-Robertson- Walker (FRW) space-time and 
all homogeneous anisotropic modifications thereof. The invariant 3-manifolds of Bianchi-type VIII are 
non-compact and describe open space-times. (The open FRW space-times are not of this particular 
Bianchi-type, however, see Ryan and Shepley 1975 |4[). 

We shall restrict the further discussion to diagonal models, where Qijif) can consistently be chosen 
diagonal, 

g = J-diag ( e (2a+2/3 ++ 2V3/3-) ; e (2 a+ 2^ + -2V3/3-) ; e (2a-4/3 + )\ (5) 
67T ^ ' 

with three time-dependent parameters a, /3 + , [3-. The special space-time metric ([[]) can now be inserted 
in the action functional of general relativity, in particular in its first-order (Hamiltonian) form first given 
by Arnowitt, Deser, and Misner 1962 [[37]. We shall neglect the matter contribution in the action, which, 



in the anisotropic case, can be shown to become unimportant close to the initial singularity. The result 
is a reduced action of the form 

J d\(p v q v -N(X)H(q,p)) (6) 



S 
with 



N{X)d\ 




? = (7) 



H(q,p) = ^(- P l+pl+ P 2 _)+V(a,f3 + ,(3„). 
Here N(X) > is a Lagrange parameter called the lapse function, whose variation enforces the condition 



H = 
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(8) 



q u , p v are used to denote the coordinates in the space of diagonal, homogeneous 3- metrics (mini- 
superspace) a, /?+,/?_ and their canonical momenta. The action (^) has the same form as for a particle 
in a 3-dimensional potential well, except for the Minkowskian nature of the metric in mini-superspace. 
The potential V(a,/3 + ,j3-) is given by the scalar curvature of the homogeneous 3- metric 

V = -12ir 2 det(g ij ) ( - 3) R. (9) 

For the two Bianchi-types of special interest here it is given by (Ryan and Shepley 1975 M) 



v = y 



2e 4 ^+ (cosh (4y/30-) - 1) + e~ 8f3 + 
±4e" 2/3 + cosh(2 v / 3/5-) 



(10) 



where the upper and lower sign apply to Bianchi-type VIII and Bianchi-type IX, respectively. 

The action (§) with eqs. flSp-fllOD correctly generates the reduced equations of motion obtained 
from Einstein's equations with the ansatz (^). It therefore gives an exact description of the spatially 
homogeneous models. 

We note that the potential flUf ) for Bianchi-type VIII is positive everywhere and decreases towards 
for /?_ = 0, (3 + — > oo. It is therefore like a space-time dependent mass of a relativistic particle with 
a playing the role of the time-coordinate changing monotoneously with A. For Bianchi-type IX, on the 
other hand, the potential is everywhere increasing for \f3 + \ —>■ oo or |/3_| — > oo, and V becomes negative 
in a neighborhood of the origin /3 + = /?_ = 0. Therefore the action of V cannot be compared to a 
mass-term in this case, and a need not change monotoneously with A. 

However, it can be shown (Bogoyavlensky and Novikov 1973 0, Rugh 1992 |l(J) that a cannot have 



a minimum fc = 0, > and can have only a single maximum. To see this we use H = (and N = 1) 
to express 

d 2 a dH . _ _ n fda\ ^2^2 

-l^ = -^-- 4V -- 2 {Tx) + + (") 

which implies d 2 a/d\ 2 < (except for p + = p_ = 0)for da/dX = 0. Hence, sufficiently close to the 
singularity a — ► —00, where a maximum of a does not occur, a or any monotoneous function of a can 
in fact be used time-coordinate. 

For a — ► —00 the potential ( |10|) approaches zero for values of (3 + , /3_ inside a triangular region of the 
-plane bounded by the three straight lines 

a + ft. ± v^ft = 
a - 2(3+ = 0. 

The three sides of this triangle recede with decreasing 'time' a with velocity 1/2. The velocity of the 
system point with respect to the time a in the region where V is negligible is approximately 1. The 
potential (|10D rises steeply (exponentially) outside the lines (0) with the exception of narrow channels 
reaching to infinity. These channels are located symmetrically at all three corners of the triangle for 
Bianchi-type IX and in them the potential approaches zero from below for a — > — 00. For Bianchi-type 
VIII there is only one channel at the corner /3_ = 0, (3 + = —a reaching to (3 + = +00, and there the 
potential approaches zero from above. In an apparently reasonable approximation in the limit a — ► — 00 
the exponentially rising potential walls can be replaced by the infinitely steep walls of a billiard located 

TV 



at the straight lines (0). The picture of a billiard ball of unit velocity elastically bouncing from the walls 
of an equilateral triangle receding with velocity 1/2 permits a thorough and lucid qualitative discussion 
of the resulting dynamics, which has been given in many papers and textbooks, starting with the classic 
works of Belinski et al 1969 |IJ and Misner 1969 ||. Detailed numerical simulations fully confirming the 
qualitative picture have recently been reported e.g. by Rugh 1992 |TTJ. We refer to the quoted literature 
for details. The channels at the corners, and the (for a — > — oo increasingly) unlikely events, where 
the system point enters one of those channels, are not described by the billiard approximation. This 
fact makes it hard to assess the quality of the billiard approximation. In particular, chaos in the billiard 
approximation need not imply chaos in the unapproximated system, whose presence or absence is, in fact, 
still an open question (see Cantopoulos et al 1993 |]TTJ for an indication that the system might be actually 
integrable; however, it should be noted that eq. ([□]) e.g. implies that the phase-space function p a e~ 2a 
varies monotonously with A which makes the system very special indeed, and calls for extreme caution 
before drawing too strong conclusions from the positive result of the Painleve test for integrability) . 

In the next two sections we shall be concerned with the properties of the cosmological models in 
the billiard-approximation. However, in the final section we shall return to the description by the full 
Hamiltonian (8). The 'free'motion between two bounces is called a 'Kasner epoch' because the Kasner 
universe is described by the Hamiltonian (8) with V = (Ryan and Shepley 1975 [§]]). Repeated bounces 
between the same two walls constitute an 'era'. For a — > — oo , the motion then consists of an infinite 
sequence of 'eras' following each other in which the two walls defining each era are permutated in a 
chaotic manner (Barrow 1982 [|, Khalatnikov et al. 1985 0). 



REDUCTION TO A BILLIARD ON THE PSEUDO-SPHERE 

Discussion of chaos in the triangular billiard with moving walls are hampered by the fact that, with the 
use of a as an effective time-parameter, the system is explicitely time-dependent and the long-time limit 
a — > — oo, depending on its precise definition, is either trivial (free particle) or it does not exist (the 
particle-bounces never disappear). In fact an obvious but important point emphasized by Rugh 1992 [|K| 
and Pullin 1990 is the proper choice of the time-coordinate in discussions of chaos in general relativity: 
as Lyapunov exponents and Kolmogorov-Sinai entropy are not invariant under the free redefinitions of 
the time-coordinate permitted by general relativity some suitable coordinate condition must be posed in 
addition in order to preserve the usefulness of these mathematical concepts. Obviously one useful choice 
is a time variable which eliminates the explicit time-dependence at least asymptotically for sufficiently 
long time. 

For the problem at hand such coordinates were provided by Misner 1972 [fl^| , and first put to use in 



Chitre's thesis 1972 ]y|. These coordinates (£,£,0) are defined by 

/?+ = e* sinh £ cos 4> 

(3. = e* sinh £ sin (13) 
a = — e*cosh£ (a < 0). 

Here t is not to be confused with the synchronous time-coordinate in eq.(|l]). Rather, in the Minkowskian 

metric of the Hamiltonian (8) it is a time-like coordinate like a. In fig. 1 we give a schematic plot 

illustrating the coordinate change. The triangular billiard walls (|I2|) are shown in the (/?+, /3_)-plane for 

fixed a < 0. The corners of the triangle are on the circle f}\ + (5 2 _ = a 2 . A fixed value of £ corresponds to 

v 



Fig. 1: Billiard walls (12) and a circle £ = const in the (/?+, /3_)-plane for fixed a. 



the smaller circle 0L = a 2 tanh 2 £. The corners of the triangle are therefore at |£| = oo in the new 
coordinates. The narrow potential channels at the corner of the Bianchi VIII and at all three corners of 
the Bianchi IX potentials can therefore not even be described in the new coordinates. 
The billiard walls ([12]) in the new coordinates are at 

tanh£ = — - sec ^</> + m— (14) 

with m — 0, ±1. They are independent of the new variable t which will be used as a new time coordinate 
instead of a. The Hamiltonian expressed in the new coordinates and their canonical momenta takes 
the form 

H = \ {-V\ + V\ + - J^) e" 2 * + V& 0, t).) (15) 
Replacing the lapse function iV(A) by 

-AT (A) -> iV(A)e 2 * (16) 

the Hamiltonian becomes 

H = \ [~v\ + Pi + ^) + e 2 V(e, 0, t). (17) 

It still has to satisfy the condition H = 0. For t — ► oo the potential approaches like exp(— const. e t + 2t), 
const. > 0, inside the billiard walls (|14"D, and it approaches infinity outside these walls. Therefore 
the Hamiltonian ([T7|) ceases to be time-dependent in that limit and the canonical momentum p t is 
asymptotically conserved. The motion in this limit therefore becomes that of a free particle moving with 
arbitrary conserved energy |p 2 on a 2-dimensional space with metric 

dl 2 = rf^ + sinh 2 ^ 2 (18) 

and bouncing elastically from the infinitely steep walls at the lines (0). The Riemannian space defined 
by eq. flUf) has constant negative curvature, with Gaussian curvature K — — 1 and Riemannian scalar 
curvature R = —2. Replacing £ by r = tanh£/2 it is represented by the pseudo-sphere r < 1 with 

ci/ 2 = dr2 + r2 f (19) 
1 — H 
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or by the Poincare half-plane Imz > with 

y/3 i _ re ^+|) i ^ 

rfZ 2 = ^ 2+ 2 ^ . (20) 

r 

It is obvious from eqs. (|19D, (EOT) that the metric in these coordinates is conformal to a Euclidean metric, 



i.e. all angles are as in a Euclidean space. The walls (14) are geodesies of the metric and define a 
triangular billiard with corners at infinity. Therefore the billiard is non-compact. However, its total 
metric area is finite. In the coordinates ( [2Cj) the walls are at 

x = -l,x = 0,y>Q;(x+^ +y 2 = ^,y>0. (21) 

In these coordinates geodesies are semi-circles in the upper half of the z-plane with centers on the real 
axis. Due to the metric points on the real axis and points with y = oo are at infinity. For many 
purposes a geodesic which is repeatedly reflected on the walls is more conveniently represented on a 
somewhat enlarged domain, again bounded by geodesic walls, which are identified in opposite pairs 
corresponding to periodic boundary conditions. The smallest such domain preserving the geodesic is 
called the fundamental domain. It is shown in fig. 2 where the original billiard is displayed together with 
its extension to a fundamental domain. The letters a and b indicate which sides are to be identified. 
The conformal transformations identifying opposite sides together with the unit transformation generate 
elements of the 'Fuchsian group' of the billiard. Explicitely, these transformations are 

A : z' = z + 2; B : z' = (22) 

2z + 1 K J 

and their inverse. In order to return from the fundamental domain of the Fuchsian group to the original 
domain of the billiard, points connected by a reflection on the wall x = 0(z' = —z*) are identified. 



Chitre 1972 [TJJ first used this construction to demonstrate ergodicity and mixing of the cosmological 



billiard invoking theorems by Hopf 1936 flnj and Hedlund 1939 [B7| on geodesic flows on the quotient 



space of the pseudo-sphere and a Fuchsian group. For a similar non-compact billiard whose fundamental 
domain is the modular domain \z\ > 1, — | < x < |, y > with A : z' — z + 1, B : z' — — 1/z Artin 1924 
[ 15] already gave a demonstration of ergodicity. He used the theory of continued fractions to show that the 
transformations of the Fuchsian group, applied to the real footpoints of a geodesic semi-circle, generate 
a sequence of footpoints which lie dense on the real axis (and hence generate semi-circles which lie dense 
in the modular domain) except for a set of footpoints of measure 0. The transformations ( |22|) restricted 
to the real axis play an equally important role in the quantized billiard (see below). 

Due to the fact that the walls are geodesies the Lyapunov exponent A and the Kolmogorov-Sinai 
entropy h immediately follow from the equations for the geodesic deviations and are given by A = h — \pt\. 



QUANTIZATION 

The canonical quantization of the Hamiltonian (7) with the condition ([8]) yields the Wheeler DeWitt 
equation 

IK^ 2 - ' J- )2 - < J:' 2 ) + v(a ' W.w = o- (23) 
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Fig. 2: The cosmological triangular billiard (full lines) in the complex z-plane and its extension to a fundamental 
domain (dashed lines) with sides a, a' = A(a) and b,b' = B(b) pairwise identified. 



Eq. (|23|) is particularly simple because the lapse function was chosen so as to make the metric in the 
variables (a, /3+,/3_) flat. For a different choice of the lapse function the metric is no longer flat and 



its scalar curvature should appear in fl23]) together with the Laplace-Beltrami operator to render (|23|) 
invariant under such conformal transformations (Misner 1972 fl2f ). In the present section we wish to 
consider eq. fl23|) in the billiard approximation where it reduces to the quantized version of eq. (IT 



(whose lapse function is still the same as in (p3|) ) with infinitely steep potential walls, i.e. 



-it 











1 d f 
sinhf 1 d£ 



d 2 



ip = o 



(24) 



dt V dtj sinh f d£ \~~ ~ ' 0^ J sinh^ ^ , 
with Dirichlet boundary conditions on the lines ([14]) . Separating the coordinate t from ^, by the ansatz 



^(t,e,0) = e" ( ^ + ^V(e,0) 



and introducing the coordinates ( p0|) we obtain 



2 



+ 



5 2 



(25) 



(26) 



with Dirichlet boundary conditions on the lines (|2l|) . We note that the eigenvalues \u\ give the allowed 
values of the observable \pt\ which classically equals the Lyapunov exponent and the Kolmogorov Sinai 
entropy. 

There is an immense mathematical literature on the eigenvalue problem (p6|), see Hejhal 1976 [39], 
Helgason 1981 |40| , Venkov 1982 [II]]. One important result is that the eigenvalues u of eq. ( ^6|) are 
given exactly by the zeros s = \ + iu of Selberg's zeta function, which is defined as an infinite product 
over the primitive classical periodic orbits, labelled by p, of geometric length l p 



Z{s) 



nna 

P k 



e k (p) exp(-(s + k)l p )) 



(27) 
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where ek(p) counts the number of reflections on the Dirichlet boundaries ( pip and is given by ek(p) = 1 
if that number is even, and €k{p) — ~{~^) k if that number is odd. 

Another interesting observation, made by Bogomolny and Carioli 1992 ||28|| , is the fact that the 
solution of eq. ( |2"6| ) can be reduced to the eigenvalue problem of a 1-dimensional quantum map. In 
general, the solutions ip(x,y) of eq. have a representation in terms of functions f(t) on the real line 
(Helgason 1981 gg) 

^(x,y)= f°° f(x + t)dt. (28) 



which, like the ip, transform according to an irreducible representation of the Fuchsian group T: ip s (G{z)) = 
x(G)i/; s (z); f(G(i))\DG(t)\ 1 ~ s = x(G)f(t); GeT. Demanding, e.g., that ip transforms according to the 
unit representation of the Fuchsian group we find from (|22|) 



f(t) = /(t + 2), /(*) = /(* -2) 

/(*) = f(-^Ar)\^ + M~ 1+aiu ( 29 ) 



f(t) = /(^)|l-2t|- 1+Ww . 
The Dirichlet boundary conditions for i[) are satisfied if f(t) is restricted to odd functions 

/(*) = -/(*) (30) 
Thus, uj and fit) are determined as eigenpair from the equation 

/(*) = / (^l( mod2 )) \ 2t + l \~ X+V " ( 31 ) 

with the boundary condition (f30l). The 1-d map t' = t/(2t+ 1) (mod 2) appearing in eq. ( f3l|) is similar to 
the Gauss map. It appears as a one-dimensional map of the classical dynamics if the transformations A, B 
eq. (p2|) are restricted to the real line. The 'eigenfunctions' f(t) are rather singular objects (Helgason 1981 
flop which makes it preferable to study the adjoint map induced by eq. fl3lD on the adjoints of fit) in a 
suitable scalar product. I refer to the papers by Bogomolny and Carioli 1992 p9| , Mayer 1991 f42fl , and 
Eisele and Mayer 1993 [0 for further discussion. 



The methods described so far are interesting in their own right, but they are not the most useful ones 
if one wishes to determine a large number of eigenvalues with rather high precision. This is done best by 
solving eq. (]26f) in a suitable representation directly on a computer (Csordas et al 1991 f27H , Schmit 1991 



|3lH). The general solution of eq. fl26|) satisfying Dirichlet boundary conditions on x = —1,0 may be 



represented as 

oo 

ifrufa y) = Yj A n sm('Knx)y/yK iU} (7rny). (32) 



n=l 

Here 



2 

The Dirichlet condition on the remaining side {x + h) 2 + y 2 = \ yields 



oo 

KUz) = ^ / e' zcosht+lujt dt. (33) 



.4 2 
n=l » ^ z 




£ A n sm{<Knx)$ -- (x + l -YK luJ I - (x + i) 2 ] . (34) 
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The right-hand side can be expanded in the set {sin(7rnx)}^ =1 which is complete in [—1, 0] and we obtain 

oo 

A = E hnK = (35) 

with 



n=l 



f° 4/11 

l in = 4FiG n J ^ dx sm(nlx) sin(7ma;) W - - (x + -j) 2 K k 



{7rn\ 



4 



[x+ 2 )2) 



where Fi, G n can be chosen arbitrarily, not zero, to simplify the eigenvalue condition 



detl 



In 



0. 



(36) 



(37) 



For the evaluation of the matrix elements fl36|) asymptotic expansions of the Hankel functions are used, 



which are given in Csordas et al 1991 [27]. The eigenfunctions of the Hamiltonian transform according 
to the irreducible representations of the triangle group C% v , which has two 1-dimensional irreducible 
representations corresponding to non-degenerate states and one two-dimensional representation. For the 
purposes of level statistics each symmetry class must be considered separately. Eigenfunctions transform- 
ing under the unit representation satisfy Neumann boundary conditions on the symmetry axes, those 
transforming under the second one-dimensional representation satisfy Dirichlet conditions there, instead. 
For these cases the original triangle (|2T| ) is replaced by the smaller triangle 



x 



X 



0, y > 0; x 2 



y 



< x < o 



i.e. in the Fourier expansions 



T ->*-->-■* > 2 
-(B3p n, I on the right hand side is replaced by 2n, 21 and 



(3* 



[x 



(39) 



In this way we restrict ourselves to a desymmetrized component of the original energy spectrum. 

Remarkably, the only reason for the non-integrability of the Schrodinger eq. (p6[) , (|3~8|) is the form of 
the boundary at y — yl — x 2 which breaks the translational invariance in x-direction. It is interesting 
therefore to unfold the billiard problem and to consider not only the boundary (|39D but also the shifted 
boundary 



V = 9(cp,x) = 1 h 



\ 



\ - x 2 
c 



(40) 



depending on the parameter c p , < c p < 1. The boundaries interpolate smoothly between the 'integrable' 
boundary y = 1 for c p = 0, which approximates (38) in the interval — \ < x < reasonably well but 
preserves translational invariance, and the real boundary fl38D for c p — 1. The integrable case c p = has 
been investigated in Graham et al 1991 [26~|, the general case c p < 1 and c p = 1 in Csordas et al 1991 



27] and, with emphasis on the transitions between the 

different 'universality classes', further in Csordas et al 1993 |43 |. We shall now proceed to discuss the 



level statistics obtained in Csordas et al 1991 27 



LEVEL STATISTICS 



A very useful check on the completeness of the energy levels A = u 2 + | obtained from eq. Q57 ) up to a 
given energy A = \i is provided by the smoothed level density as obtained from the semi-classical Weyl 
formula and its refinements. For the case at hand the Weyl formula was derived in Csordas et al 1991 
[27]. If N(fi) is the smoothed number of eigenvalues A with A < fi, the Weyl formula reads 



4vriV(/i) = Afx- y/JI\n fi - By/JI + C + 0(fx~ 1/2 ) (41) 

with the area A and the constants B, C given as explicit functions of c p . The term — ^//Jln \i is due to the 
diverging circumference of the billiard. The numerically 'measured' function N(fi) deviates from N(fi), 

NQj.) = NQi) + ±-5N0m). (42) 

In figs. 3 a,b, we plot SN(fj,) for c p — 1 and c p — | for about 1400 levels, respectively. It can be seen that 
SN(fj) fluctuates around 0, indicating that no level has been missed (a systematic shift of the average of 
5N(fi) by 1 would be very noticeable in fig. 3). Another immediate observation from fig. 3 is that the 
fluctuations are significantly larger in the 1 than in the CclSG Cp — 1/2. Using the smoothed 



Fig. 3: The deviation of the spectral staircase from the Weyl formula (a) c p = 1 and (b) c p = 
et al H3). 
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i (from Csordas 



level-staircase (£T) the average level density can be scaled out by the transformation 



Xi = N(\i). (43) 

Then it is possible to determine the nearest neighbor level-spacing distribution P(s)ds which gives the 
probability to find the next level Xi + \ in the interval (xi + s, X{ + s + ds) if there is a level at Xi. The 
integrated level spacing distributions I(s) = J S P(s)ds for plotted in fig. 4 a,b, respectively, 

and compared with the Poissonian result (dashed curve) 

I P (s) = 1 - e ~ s (44) 

and the Wigner surmise for the Gaussian Orthogonal Ensemble (GOE) of 
random matrix theory (dashed-dotted-curve) 

l w (a) = 1 - e— 2/4 . (45) 

For c p = 1/2 there is good agreement with the GOE result, indicating that the transition from the 
integrable to the chaotic case has occurred. However, the case c p — 1 fits the Poissonian distribution 



better than the GOE distribution, which was surprising at first (see Balazs and Voros 1986 [17| for the 
first observation of this phenomenon in a compact billiard), indicating some hidden non-generic behavior 
in this prototypical chaotic system. This non-genericity is now understood to be a consequence of hidden 
degeneracies occurring in some billiards on the pseudo-sphere which are due to certain special arithmetical 
properties of the underlying Fuchsian group (see below). 
In fig. 5 a,b we also plot the averaged spectral rigidity 



1 r x+p/2 

A 3 (L,x) = - A 3 (L,x')dx' (46) 

P Jx-p/2 



with 

1 rx+L/2 

A 3 (L, x) = min - / (N(x) - Ax' - Bfdx' (47) 

A,B L Jx-L/2 

again for c p = 1,1/2. For comparison the Poisson statistics A3 = L/15 and the GOE statistics A3 ~ 
\ InL — 0.00695 . . . for large L are also shown. For c p — 1/2 the GOE result fits the data well, while for 
c p = 1 the As-statistics for small L is again closer to the Poissonian result. For larger L the Poissonian 
result is not expected to be a good approximation, not even in the integrable case, because of the influence 
of periodic orbits (Berry 1985 [24]). For the integrable case c p = see Graham et al 1991 PEfl . The 



conclusion that the case c p = 1 is non-generic and similar to an integrable case is further confirmed by 
the S 2 -statistics 

S 2 = (K(r) - (n t (r))} 2 ) (48) 

where Hi is the number of levels in an interval of size r. These data are shown in figs. 5 a, b respectively. 
An explanation for the strong deviation from GOE statistics in the chaotic case c p = 1 has been given by 
Bolte et al 1992 [3T1] and Bogomolny et al 1992 [p8[| . It is due to the existence of further symmetries of 



some of the Fuchsian groups (called 'arithmetic'), to which the infinite equilateral triangle also belongs. 
Classically these additional symmetries give rise to the appearance of a multiplicity (g) of periodic orbits 
of the same length / with (g) for length I < M growing exponentially with M/2 for large M 

(9) = = const ^rr- I 49 ) 
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When used, in Berry's 1985 semi-classical theory of level statistics the large value of (g) for M » 
1 has consequences very similar to those of the infinite multiplicities of periodic orbits in integrable 
systems. Quantum mechanically for arithmetic Fuchsian groups there exists an infinite denumerable 
set of hermitian operators T n 'Hecke operators', commuting among themselves and with H, which can 
be diagonalized simultaneously with H (Hejhal 1983 J5D , Venkov 1990 Again one then expects 

Poissonian statistics, as observed. Even a slight shift of one of the walls (c p < 1) is enough to destroy 
the length-degeneracy of the periodic orbits (or the conservation of the Hecke operators) and leads to a 
transition to GOE statistics for sufficiently high-lying levels. This transition has recently been studied 
in detail by Csordas et al 1993 [Q, who found a scaling behavior in the variables 1 — c p and the energy 
quantum number n. 



Fig. 4: The integrated nearest-neighbor level spacing distribution from approximately 1400 levels for (a) c p = 1 
and (b) c v = \. Solid, dashed, and dotted lines are the numerical result, the Poissonian, and the GOE predictions, 
respectively (from Csordas et al p7[). 
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Fig. 5: The spectral rigidity A^(L,x) at x = 700 for (a) c p = 1 and (b) c p = \. Solid, dashed, and dotted lines 
are the numerical result, the Poissonian, and the GOE predictions, respectively (from Csordas et al p7||). 



SEMI-CLASSICAL SOLUTIONS OF THE WHEELER DEWITT EQUATION 



The study of the separated eigenvalue euqation (|26|) presented in the two preceding sections has several 



drawbacks: (1) It can claim validity only in the asymptotic region t — > +oo where t is the function 
of a,j3l_ + fi 2 _ defined in eq. (13). (2) It completely neglects the influence of the potential channels at 
the corners of the triangle in configuration space, which become extremely narrow for t — > +oo, but, in 
principle remain always open. (3) Even if eq. (|26|) gave a complete description it would leave open the 
question which linear superposition of the eigenstates of (^6|) determines the actual wave-function of the 
Universe. 

For all these reasons it seems desirable to avoid the transformation (13) and instead to go back to 



the original Wheeler DeWitt equation (|23"D , even at the cost of treating a more difficult problem. 

As we are still interested in the asymptotic region where the influence of matter is negligible we may 
anticipate that the solution of (p3|) will have the form of a bound state 



1>(a,P + ,l3-)~exp(-<l>(a,!3 + ,l3-)/K) (50) 
where we restore Planck's constant here and in (|2"3"D in order to have a formal control over the semiclassical 



limit h — > 0. Surprisingly, exact analytical solutions of the Wheeler DeWitt equation of the form (^) can 
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Fig. 6: £2(0 statistics for (a) c p = 1 and (b) c v = \- Solid, dashed, and dotted lines are the numerical result, 
the Poissonian, and the GOE predictions, respectively (from Csordas et al [27]). 



be given (Graham 1991 |54|, Moncrief and Ryan 1991 p3| Bene and Graham 1993 fl35|). The potential 
in eq. ( |50"D may be interpreted as a Euclidean action which is required for the system in its bound state 
to tunnel from some initial point a , /3 +0 , /3-o to a given end-point a, /?+,/?_. (In this context it may be 
helpful to recall a harmonic oscillator with coordinate x in its ground state where in appropriate units 
= x 2 is the Euclidean action for tunnelling from x = to the classically forbidden amplitude i^O). 

Inserting the ansatz fl50|) in eq. (p3|) (and taking ft — > or adding suitable terms of order ft to 
the potential V which is always possible by an appropriate choice of factor ordering when defining the 
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quantized model, See Moncrief and Ryan 1990 |35|) we find the Euclidean Hamilton Jacobi equation [] 

2 - / n i \ 2 , / n , \ 2 

-V(a,/3+,/3_) = (51) 

which now contains the inverted potential —V. The following exact analytical solutions of this Euclidean 
Hamilton Jacobi equation are known, Belinsky et al 1978 |[46|| , Gibbons and Pope 1979 ^7|, Atiyah and 
Hitchin 1985 | i8| , in the case of Bianchi type IX 



1 / g 2/3 + +2v / 3/3_ + e 2/3+- v / 3/3_ + e ~4f3+ 




_ e -^ ++ ^_ _ e -^ + -V3/3_ j _ (52) 

From 03 two other solutions may be obtained by 120° degree rotations in the (j3 + , /5_)-plane. It is clear 
that — is also a solution if is a solution. In the following we shall only consider the solutions 0i,02 
which have the full triangular symmetry of the problem. In figs. 6 and 7 we give plots of the wave- 
functions ( |50"1) obtained with 0i, 02 for different fixed values of a. The maximum -0 is always normalized 
to unity. It is seen that the channels indeed modify the wave-function. However, because the <pi are real 
there is only a tunnelling penetration into the channels without any real propagation inside them. 

For the case of Bianchi type VIII very similar analytical solutions exist. They are given by (Bene and 
Graham 1993 g|) 

1 = I e 2" L2P++2V3P- _|_ e 2/3+-2 v / 3/3- _ e ~4,P+\ 

02 = ^e 2a (Ae w+ sinh 2 V3(3_ - e~ 4/3 + ± 4ie~P + cosh V3(3_) (53) 

03 = ^e 2a (4e 2/3+ cosh 2 ^3/3- - e" 4/3 + ± 4ie" /3+ sinh V3/?_) 

The fact that 02 and 03 turn out to be complex is due to the non-binding nature of the Bianchi VIII 
potential which now makes a real propagation of the wave-function (without tunnelling) inside the channel 
along the positive /3 + -axis possible and gives rise to a real physical action (i.e. an imaginary Euclidean 
action). For a — > — oo and values of inside the triangle (|II]) the imaginary part of is negligibly 

small. Yet solutions of the form ( |50"D with = ±0^, 0j of eq. (53), turn all out to be non-normalizable 
for fixed a, i.e. none of these solutions is physically acceptable for Bianchi type VIII. 

We now turn to the problem of physically interpreting the special solutions (|52"D- First of all, this 
involves the question of the interpretation of the wave-function of the Universe, which would need a 
discussion we don't wish to enter here (see e.g. the discussion in Bene and Graham 1993 [PH]). We 



*A Euclidean Hamilton Jacobi equation (|5l| ) with a potential Vo(a, /3+, /?_) via eq. ( pty) corresponds to a Wheeler 
DeWitt equation (p3|) with a potential 



V = V (a,/3 + ,P-) + 



h fd 2 4> d 2 cp d 2 



2 V<9a 2 d(3 2 _ d(3 2 _ 



Such a potential appears naturally in supersymmetric formulations (Graham 1991 [p3[ , Bene and Graham 1993 

However, in any case V is reduced to Vo in the semiclassical limit 7i — > 0. 

YVT 



Fig. 7: Wormhole state: The wave function |^| 2 for <\> = <\>\ as a function of /3+,/3_ for four different values 
of a. The maximum of |0| 2 is always normalized to 1. The changing scale of /?+,/?_ should be noted (after 
Graham 1991 @). 
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Fig. 8: No-boundary (Hartle- Hawking) state: The same as Fig. 7 but for <fi = 4>2- 
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shall simply interpret \ip\ 2 for fixed logarithm of the scale parameter a, as the probability distribution 
of the anisotropy parameters P+, P- for that scale parameter. However, we still have to understand 
which wave-function applies to which physical situation. This question can be answered by recalling the 
meaning of as a Euclidean action. The (imaginary-time) solutions of the classical equations of motion 
corresponding to these actions can be interpreted (for h — > 0) as the most probable paths followed in 
the tunnelling propcess. These paths, in the present case, correspond to a sequence of 3-spaces followed 
in the tunnelling process from a given intial 3-space up to the final 3-space whose parameters a, 
appear as the arguments of 0. The initial point of the tunnelling path is chosen as a most probable state 
=min for a fixed where the system likes to sit when it does not fluctuate away on one of its tunnelling 
excursions to the neighborhood (again, recall the harmonic oscillator in its ground state for a simple 
example of this physical picture). Let us now consider the most likely tunnelling paths. They satisfy 





da 


Pa 


~dX = 




d/3 + 


P+ 


dX 


P- 


d(5_ 


dX 


d(3 + 


1 d 


da 


2 dp, 


dp_ 


1 d 


da 


2 dp. 



da 



(54) 



dp. 

It follows that 

* -u(p+,P-) 

■U(P+,P-) (55) 

aa z cp_ 

with the new a-independent potential 

U{P + ,P_) = \n\e- 2a <P{a,P + .p.)\. (56) 

According to eq. (]55|) the evolution of P+,P- with increasing a is on lines of steepest descent of the 
potential U. 

Let us now consider the concrete Bianchi type IX solutions 0i and 02 of eq. (|52|) . For = 0i eqs. 
(|55l) can be solved explicitely (Belinsky et al 1978 f46fl ) yielding the Riemannian space-time metric 

,2 dp 2 1 6 /T^7~T / u 1 ^ 1 ^ 2 ^ 2 ^ 3 ^ 3 \ f „„s 

ds -^ + -/^{^ + 7^4 + 7^i) (57) 



with 

where oi, a^, are constants of integration which must be chosen in such a way that the spatial 3- metric 
of eq. fl5"7j ) becomes equal to the Bianchi type IX 3-metric with given a,P+,P- for a suitable choice of 
p. It follows from fl5"7|) that the Euclidean coordinate time (it) is d(it) = (F(p))~ 1 ^dp. The solution 
(|5"7|) must be restricted to p 2 > max(a^, a 2 , a 2 ) otherwise the 3-geometry is separated from the initial 
3-geometry by a singularity. 

The tunnelling path starts at the minimum of 0i for fixed a, which is at P + = = corresponding 
to p = oo in eq. (pTj). The given Bianchi-IX metric corresponds to the inner boundary at small p 2 (but 
p 2 > max(a^, a 2 , af)). For p — > oo the 4-metric (|57|) approaches ds 2 = dp 2 + ^p 2 (cu 1 uj 1 +uj 2 uj 2 +uj 3 uj 3 ),i.G- it 
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becomes Euclidean and flat (Belinsky et al 1978 indicating that the wave-function under discussion 
here describes the spontaneous quantum fluctuation of a flat 3-metric to one with an anisoptropic inner 
boundary or throat, i.e. a quantum wormhole (D'Eath 1993 [3D], Bene and Graham 1993 



Turning now to = 02 we note that a transparent closed-form solution like fl57p is not available in 
this case. However, we may use eqs. and the first of eqs. (54) to obtain the complete physical 

picture. Let us note that for = 2 the equi-potential lines of 02 and of U(/3 + , (3-) coincide with those of 
V(a, /3+, f3J) for fixed a. The minimum of 02 for fixed a is again at /?+ = = (3- defining the most likely 
starting point of the tunnelling path. This point corresponds to a local isotropic maximum of U, i.e. if 
/3+,/3- beginn to tunnel to small non-zero values, a increases from its inital value — oo meaning that the 
Universe expands from nothing in some random anisotropic way. From the first of eqs. (54) it follows 
that this expansion occurs with increasing A, since 02(a,O,O) < 0. However, ?7(/3+, /?_) has a degenerate 
minimum on the line where V(a, f3 + , f3J) (and hence also 02) change sign. In fact this line is formed by a 
triangle of three hyperbolas which are asymptotic to the positive f3 + axis and the two axes obtained by 
±120° rotations. Outside this triangle of hyperbolas the equipotential lines of V, 02 and U are no longer 



closed due to the existence of the channels on the three axes. Linearizing eqs. (pSD in the neighborhood 
of the line of minimal U one finds that with increasing A the line is in fact crossed, but while crossing 
a reaches a maximum and decreases again as P+,/3- start climbing the potential outside the 

triangle of hyperbolas to their given final value. 

Thus for given values of (f3 + , /?_) outside (or inside) the said triangle of hyperbolas the most probable 
path of tunnelling starts with vanishing 3-volume (a = — oo) and goes (or, respectively, does not go) 
through a maximum of the 3-volume until it reaches the given values «,/?+,/?_. The given 3-metric in 
this case forms the outer (and in fact the only) boundary of the compact Riemannian space-time region 
traced out by the 3-geometries along the most probable tunnelling path, and 02 is the Euclidean action 
of this compact space-time region. Hartle and Hawking 1983 |T{| made the general proposal that in the 



semi-classical limit the wave-function of the universe is given as the exponential of the Euclidean action 
constructed in this way. Thus 

if) = const exp(— 2 (a, /3 + , (3J)/K) (59) 
is identified as the no-boundary state of the Bianchi-type IX quantum space-time (Bene and Graham 1993 



|35|). The state has also been discussed by Hawking and Luttrell 1984 [TZOJ and by Moss and Wright 1985 
however, this work was qualitative or numerical and the analytical expression (|59| ) was not given 
there. 

A similar discussion for 3 of eq. ([52]) shows that the wave-function (50) for this case arises by quantum 



tunnelling from a disk shaped 3-geometry (/3 + = oo, (3_ = 0) with infinite volume (a — > oo) down to 
a given 3-geometry a, /3 + ,j3- with a decreasing monotonously along the most probable tunnelling path 
which closely follows the /3 + -axis before it finally turns to the given value /3_ . This is again a solution of 
wormhole type, but this time of a virtual wormhole in a highly anisotropic 3-geometry. 

Finally, let us compare the solutions obtained in the present section with the solutions obtained in the 
two preceding sections. First, we note that the solutions given there have only restricted validity, as they 
were based on an approximation which is only valid in the limit a — > — oo. On the other hand, in that 
limit, the Wheeler DeWitt equation, has been solved there, in principle, as a general linear superposition 
of the solutions of the eigenvalue problem of the two-dimensional billiard. In order to select the special 
cosmologically relevant solution an additional 'initial' or boundary condition would still be needed. The 
simplest choice is to select the ground state of the billiard (Graham and Szepfalusy 1990 but other 
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choices may, of course, be possible. By contrast, in the present section we have obtained some particular 
solutions only, which are, however, valid for arbitrary a. As discussed above, these solutions correspond 
to the special boundary conditions of the 'no-boundary state' and the wormhole state. Their cosmological 
significance is therefore clear. If specialized to the regime a — > — oo these particular solutions may be 
represented as linear superpositions of the eigenmodes of the two-dimensional billiard, i.e. they can be 
used to determine the probability amplitudes of the eigenmodes of the billiard in the wave-function of the 
Universe. As can be seen from figs. 6, 7 for a — > — oo not only the ground state of the billiard, but also 
arbitrarily high lying states will be excited, due to the increasingly sharp drop of the wave-function from 
a nearly constant value inside the billiard to at the boundary. In fact it appears that these solutions 
asymptotically approach the constant solution of the billiard for Neumann boundary conditions which 
contains arbitrarily high lying states satisfying Dirichlet boundary conditions. The kinetic energy of the 
billiard is the positive energy residing in the anisotropy degrees of freedon /3+,/3- or x,y. This positive 
energy is exactly balanced, because of if = 0, by the negative energy — \p\ of gravitational attraction. 
The particular solutions discussed in the present section, in the limit a — > — oo, contain components for 
which both energies are unbounded. 
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